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Spatial Kerr solitons, typically associated with the standard paraxial nonlinear Schrodinger equa- 
tion, are shown to exist to all nonparaxial orders, as exact solutions of Maxwell's equations in the 
presence of vectorial Kerr effect. More precisely, we prove the existence of azimuthally polarized, 
spatial, dark soliton solutions of Maxwell's equations, while exact linearly polarized (2+l)-D solitons 
do not exist. Our ab initio approach predicts the existence of dark solitons up to an upper value 
of the maximum field amplitude, corresponding to a minimum soliton width of about one fourth of 
the wavelength. 

PACS numbers: 42.64.-K, 42.65.Tg, 42.81. Dp 



Bright and dark optical spatial solitons have been and 
still are the object of an intense theoretical and experi- 
mental investigation. 0, 0] In particular, Kerr solitons 
have always represented the theoretical reference model 
in view of the simplicity of their analytic description. 
This is mainly due to the fact that they are solution of 
an equation, the so-called nonlinear Schrodinger equa- 
tion (NLS), which is exactly integrable and also admits, 
in particular cases, of simple analytic solutions. The 
main limit of the NLS is that it is only able to de- 
scribe scalar optical propagation in the paraxial approx- 
imation. Nonparaxial descriptions have been introduced 
by adopting suitable asymptotic expansions in the small- 
ness parameter X/d (where d is the beam width) which 
result in improved versions of NLS including both scalar 
EE El and vectorial 0, H, @ contributions. Besides 
eliminating the unphysical catastrophic collapse associ- 
ated with standard NLS [HIHII3], the ^ show that the 
transverse Cartesian components of (2+l)-D propagating 
beams are mutually coupled, so that linearly polarized 
(2+l)-D Kerr solitons cannot exist. The same nonparax- 
ial approach can also be used to predict the existence of 
(1+1)-D bright and dark solitons EL to the first 
significant order in X/d. 

The standard or improved versions of the NLS are lim- 
ited by the underlying approximation scheme and fail 
whenever X/d approaches one. The question naturally 
arises: are there solitons to all nonparaxial orders? More 



in general, is it possible to find exact soliton solutions of 
Maxwell equations? 

In this paper we find, by starting from Maxwell's equa- 
tions, that an exact solution (to all nonparaxial orders) 
indeed exists in the form of an azimuthally polarized dark 
soliton, i.e. 
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where r, ip, z are cylindrical coordinates with unit vec- 
tors e r ,e v , e z and x is a real number. More precisely, we 
are able to prove the existence of azimuthally polarized 
dark solitons, i.e., of fields of the kind given in Eq.Q), 
where E(r) monothonically ranges from E(0) = to 
£■(00) = Eoa and Ex, is a given real constant. We numer- 
ically evaluate the soliton shape and its existence curve, 
our results being obtained without any approximation, so 
that the whole range of possible soliton widths is con- 
sidered without any formal distinction between paraxial 
and nonparaxial regimes. We note that the existence of 
azimuthally polarized, circularly symmetric fields of the 
form E(r, z)e v , of which the soliton given in Eq.Q is a 
specific case, is due to the simple symmetry properties of 
Kerr nonlinearity. 

A monochromatic electromagnetic field 
i?e[Eexp(— iut)], i?e[B exp(— iut)] propagating in a 
nonlinear medium obeys Maxwell's equations 
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where no labels the linear refractive index and P„i is 
the nonlinear polarizability. In the case of nonresonant 
isotropic media |l5j , the vectorial Kerr effect is described 
by the polarizability 



-eon n 2 



|E| 2 E- 



f 



-(E-E)E* 



(3) 



ri2 being the nonlinear refractive index coefficient. After 
eliminating B from Eq.J^J) we get 



V x V x E = fc 2 E + k 
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where k = uquj / c. Let us consider fields of the form 

E(r,<pz)=E< fi [r,z)e ip + E z (r,z)e z , (5) 

describing a circularly symmetric field with vanishing ra- 
dial component. Inserting Eq.JSJ) in Eq.Q, we obtain 
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Internal consistency of the set of Eqs.® (three equations 
in two unknowns) requires E z — 0. As a consequence, the 
second of Eqs.© yields 
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We note that circular symmetry and polarization im- 
posed to the field, together with the symmetry properties 
of Kerr effect, have allowed us to reduce Maxwell's equa- 
tions to the single Eq.0. Equation {7|) is conveniently 
rewritten in the dimensionless form 



where p — \/2kr, ( = kz, U = \J\n2\jnoE v and 7 
n2l\n2\- If we look for soliton solutions of the form 



where (see Eq.Q) a 
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x/k, Eq.© becomes 

l)u- 7 u 3 . (10) 
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FIG. 1: Soliton profile u(p) for various values of u a 



Both the structure of Eq. 1)1 OJI and the azimuthal field po- 
larization dictate u(0) = 0, so that azimuthally polarized 
bright solitons do not exist. In order to find dark solitons, 
we introduce the further condition 



lim u{p) = u c 

p — >OG 



(11) 



together with the vanishing of all derivatives for p — ► 00. 
Since focusing media (7 = 1, i.e., n% > 0) are not able to 
support dark solitons, we consider hereafter defocusing 
media (7 = —1, i.e., n 2 < 0), so that Ea. lfTnj) reads 



d f du u 
dp \dp p 



-(a 2 -l)u + u 3 , (12) 



which implies, together with the above boundary condi- 
tion at infinity, 



±^1-2^ 
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While positive and negative signs of a respectively refer 
to forward and backward travelling solitons (see Eq.©), 
u(p) depends on a 2 (see Eq.^HlO and is clearly the same 
in the two cases. Equation H13J) shows the existence of 
an upper threshold for the soliton asymptotic amplitude 
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since, otherwise, a would become imaginary. From an 
intuitive point of view, the existence of this threshold 
is related to the dominant defocusing effect due to the 
nonlinearity over the focusing one due to diffraction. If 
we now insert Ea. (|13fl into Ea. H12l) . we obtain 



d f du u 
dp \dp p 



( u2 - U lo) U - 



(15) 



We have carried out a numerical integration of Eq. I|15|) 
with boundary conditions u(0) = and u(oo) = Urx, 
by employing a standard shooting-relaxation method for 
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FIG. 2: Normalized asymptotic optical intensity Ioo/Io as a 
function of the asymptotic adimensional field amplitude Moo • 
Note that two solitons exist for any allowed asymptotic optical 
intensity. 



boundary value problems. The results of our simulations 
show that dark solitons can be obtained in the range 
of field amplitudes < < l/\/2. Different soliton 
profiles are reported in Fig.l. 

In order to complete our analysis, we now evaluate 
both the magnetic field and the Poynting vector. Recall- 
ing the expression of the soliton electric field 
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we obtain, from the first of Eqs.(|2J) written in cylindrical 
coordinates, 
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The magnetic field has a radial component whose shape 
coincides with that of the electric field, and a vanishing 
azimuthal component, so that E and B are mutually or- 
thogonal. With the help of Eas. lfTo)) and {T7J, the time 
averaged Poynting vector 



S = —ReCE x B*) 
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turns out to be given by 
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We note that S is parallel to the z— axis, consistently 
with the shape-invariant nature of solitons. From an an- 
alytical point of view, this corresponds to the it/2 phase 
difference between B z and E v (see Eas.lfTHf and lfTT|) ). 
As expected, the Poynting vector is either parallel or an- 
tiparallel to e z according to the sign of a, while its ampli- 
tude is proportional to |E| 2 . The above plane wave- like 



FIG. 3: Existence curve relating the normalized soliton opti- 
cal intensity HWHM to 



properties are consistent with the nondiffractive nature 
of exact solitons. 

It is worthwhile to underline that, in the case of the 
azimuthal dark solitons we are considering, the asymp- 
totic optical intensity loo = |S(oo)| turns out not to be 
proportional to u^. In fact, by using Eas. (fH5| and fTT))l . 
one obtains 



(20) 



where Iq = kriQ / {2uj po\n2\) , whose profile is reported in 
Fig. 2. Equation lj2H shows that the asymptotic opti- 
cal intensity is not a monotonically increasing function 
of the asymptotic field amplitude, but reaches its maxi- 
mum threshold value I^ ax — Io/3 3 ^ 2 in correspondence 
to = l/v3- This is connected to the a— depen- 
dence of the magnetic field (see Ea. (|17fl ) whose radial 
part tends to vanish for Uoo — » 1/V2. A related and 
relevant consequence of Ea. (|20J) is the existence of two 
solitons of different widths for a given asymptotic optical 
intensity. The existence curve relating the normalized 
half width at half maximum (HWHM) of the soliton op- 
tical intensity profile |S(p)| to Uoo is reported in Fig. 3. 
In particular, Fig. 3 shows the existence of a normalized 
minimum HWHM ~ 2.1 (~ 0.24A) for Uoc = 1/V2. In 
addition, Fig. 2 shows that a normalized HWHM ~ 2.7 
(~ 0.3A) corresponds to Uoo = 1/ \/3 for which the soliton 
attains the maximum asymptotic optical intensity I™ ax . 

It is interesting to examine the behavior of our solution 
in the limit of large p. To this end, neglecting in Eq. Ijl5(l 
the term in uj p, we have 



d 2 u , 9 
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which formally coincides with the equation describ- 
ing one dimensional linearly polarized paraxial dark 
solitons. Equation (|21|) admits of the solution u = 
Moo tanh(/9u 00 /v / 2). This solution can be compared with 
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FIG. 4: Plot of the ratio R(p) = tioo tanh(u CXJ p/\/2)/u(p) for 
different values of 

the exact one. This is done in Fig. 4 where the ratio R{p) 
between the hyperbolic tangent and the exact solution is 
reported as function of p, for different values of tioo- 

The 

hyperbolic tangent solution reproduces the exact one for 
large values of p, as expected, while it at most differs by 
a factor = 1.2 for small values of p. 

Recent findings indicate that radially and azimuthally 
polarized laser beams can be efficiently gene rated |l(i Il7| 
and focused to subwavelength diameter [lj| in resonator- 
like schemes. These, combined with an appropriate non- 
linear material can form a basis for a first experimental 
investigation of the existence of the azimuthally polarized 
solitons. As an example, if we consider a strongly defo- 
cusing nonlinear medium, like Sodium vapor at a density 
of ~ I0 12 atoms/cm 3 for which m = — 4- 10~ 14 {raj "Volt) 2 
(lif. the upper threshold for the asymptotic optical soli- 
ton intensity I^ ax is of the order of .5 MW/cm 2 . 

In summary, we have found, to the best of our knowl- 
edge, the first exact, two-dimensional, dark spatial Kerr- 
soliton. This soliton is obtained as a particular solution 
of an equation (Eq.©), which in turn exactly follows 
from Maxwell's equations, describing circularly symmet- 
ric, azimuthally polarized electric fields. Our ab initio 
approach inherently overcomes any approximated scheme 
employed to describe paraxial or slightly nonparaxial sit- 
uations. In particular, the derived analytical expression 
of the soliton propagation constant as a function of the 



normalized (p — > oo) field amplitude shows that the 
soliton does not exist whenever Uoc > l/v2- We have 
also numerically evaluated the existence curve relating 
the soliton width to Uoa] as a relevant result, this curve 
shows that the soliton width attains the minimum possi- 
ble value of about A/4 for approaching 1/^/2. 
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